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ABSTRAGT 

We present an unprecedented measurement of the disc stability and local instability scales in the 
luminous infrared Seyfert 1 host, NGG 7469, based on ALMA observations of dense gas tracers and 
with a synthesized beam of 165 x 132 pc. While we confirm that non-circular motions are not significant 
in redistributing the dense interstellar gas in this galaxy, we find compelling evidence that the dense 
gas is a suitable tracer for studying the origin of its intensely high-mass star forming ring-like structure. 

Our derived disc stability parameter Q accounts for a thick disc structure and its value falls below 
unity at the radii in which intense star formation is found. Furthermore, we derive the characteristic 
instability scale Ac and find a striking agreement between our measured scale of ~ 180 pc, and the 
typical sizes of individual complexes of young and massive star clusters seen in high-resolution images. 

Subject headings: galaxies: ISM — galaxies: kinematics and dynamics — galaxies: individual (NGG 
7469) 


1. INTRODUCTION 

The central kpc region of the barred Seyfert 1 host 
NGG 7469 displays a number of morphological and kine¬ 
matic features that provide interesting clues on the dy¬ 
namical behaviour of the interstellar medium and its con¬ 
nection with star formation across the disc in this galaxy. 


et al. 


19901): 

Salamanca et al. (1995) 

Davies et al. (2004 

2DD5) and 


ft 


icks et ai. 


plays intense star formation in ai?^ 0.5 — 1 kpc ring-like 
structure contributing to up to 2/3 of the galaxy’s bolo- 
metric luminosity (Lboi = 3 x lO**^©). The ring-like 
structure is dominated by young (< 100 Myr) and mas¬ 
sive (lO^-lO*" Mq) supe r star clusters, with the enclosed 
mass of ~ 2.7x 10® Mq {Mauder et al. 1994 Diaz-Santos 


et al. 


20071. The individual clusters are g rouped in 
150 - 2U0 p c-scale star forming complexes (Diaz-Santos 
et al. |[2007|. The relatively gas poor nuclear 30 - 60 pc 


region hosts a < 100 Myr nuclear star cluster which con¬ 
tributes to 10% of the nuclear bolometric luminosity and 
30% of the nuclear AT-band continuum light. Further¬ 
more, the intensities of the cold and warm molecular gas 
are consistent with a GO-to-H 2 conversion factor similar 
to other star forming galaxies, i.e., Xga = 0.4 - 0.8 the 
conversion factor in the Milky Way (Davies et al. ||2004 
Bolatto et al. |[^013 |. 


The dynamics of the molecular gas around the nucleus 
is consistent with a thin disc model with ~ 10® Mq 
dynamical mass and ~ 3100 Mq pc“® surface density 


(Hicks et al. 2009). The stellar structure inside the 


nuclear 30 pc is dispersion dominated whereas the Vja 
rapidly raises at this radius, staying c onsistently around 
unity out to ~ 150 pc. Furthermore, [Hicks fc Malkan~| 
(20081 applied dynamical models to measure the mass of 
the central supermassive black hole to ^ 5 x 10^ Mq. 

The molecular bar-like structure observed in GO maps 
taken with the IRAM millimeter interferometer on the 
Plateau de Bure shown no cle ar signatures of bar kine¬ 
matics (Davies et al. 2004) and it does not have a 


counterpart m high resolution V or AT-band images o r 


radio continuum maps a nalyzed byjGenze l et al.~|(|1995 ); 
Mulcha ey et al. |(|1997|);|La i et al. |(|1999D:|Scoville et al. 

[la et ai. l 200ip (see also Fig. . As- 


and Golina 

sociatmg the molecular bar-like morphology with a non- 
axisymmetric perturbation raises the problem of locating 
the position of the intensely star-forming ring-like struc¬ 
ture. With its dominant recent burst of star formation, it 


onance radius (e.g., 

Regan & Teuben |2003 

Pihol-Ferrer[ 

et al. ||2014 

Li et a 
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Given that no signilicant kinematic characteristics ot a 


it has been suggested that what appears to be a bar is 
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Fig. 1.— The central region of NGC 7469 with ellipses marking 

the 500 and 1000 pc boundaries in the plane of the galaxy disc. 


likely a blending of gas clumps. Nevertheless, these au¬ 
thors assumed epicyclic approximation to associate the 
location of the ring-like structure with the inner Lind- 
blad resonance radius of a primary bar, rotating at a 
pattern speed of < 80 km s“^ kpc“^. Accordingly, the 
pattern speed for the region inside the ring was esti¬ 
mated to 235 ± 5 km s~^ kpc~^. Congruent with this 
picture, the star-forming regions in NGC 7469 would be 
supported by its non-axisymmetric structure(s). How¬ 
ever, strong bar/spiral-induced flows would be expected 
to sustain the large amounts of molecular gas (gas-to- 
dust ratio ^ 600) with HCN J(4 —)• 3)/CO J{3 2) 

~ 1 ( Papadopoulos et al ||2012[ ). 

While these results provide important clues on the fate 
of the interstellar gas in NGC 7469, a concise picture for 
the dynamical effects that trigger the intense star forma¬ 
tion in this galaxy is missing. To build a realistic scenario 
for the relative role of the local and global instabilities 
that ignite these intense bursts of star formation in the 
central kpc radius of NGC 7469, we present here a de¬ 
tailed analysis of the distribution and kinematics of the 
densest phase of the molecular gas in the region of inter¬ 
est. We use the HCO+ J(4 —>• 3) and HCN J(4 —)• 3) 
emission lines at high spatial and velocity resolution to 
derive the instability parameter and the characteristic 
wavenumbers that describe the local instability scales. 

Throughout this work, we assume NGC 7469 to be at 
an inclination of 45° and a distance of 68 Mpc (de Vau- 
couleurs et al. |1991D , implying 330 pc/", and consistent 


with the above mentioned studies. The high angular res¬ 
olution of Atacama Large Millimeter/sub-millimeter Ar¬ 
ray (ALMA) allows us to characterize the stability of the 
molecular gas throughout the central kpc region in this 
prototype luminous infrared galaxy. We show for the first 
time, that it is possible to directly quantify the stability 
criteria that lead to the star formation in the ring-like 
structure, accounting for up to 2/3 of the galaxy’s bolo- 
metric luminosity. 


2. OBSERVATIONS, MAPS, AND PROFILES 

We observed NGC 7469 with ALMA on November 3rd 
and 4th, 2013, using 28 antennas in C32-2 configuration 
(Proposal ID 2012.1.00165.S, PI: Takuma Izumi). With 
18" primary beam size and O'.'b x Of'd synthesized beam 
size, sampled at O'.'! pixels”^, ALMA Band-7 receiver 
was tuned to cover the redshifted lines of HCN J(4 —>■ 3) 
= 354.505 GHz) and HCO+ J(4 ^ 3) (i/rest = 
356.734 GHz). The original velocity resolution was 0.43 
km s“^ but to achieve the amplitude-over-noise level > 
2, we binne d the data at an eff ective resolution of 20.2 
km s“^ (see Izumi et al. 


2015a 


for further details). 


To derive the gas distribution and kinematic informa¬ 
tion, we applied single Gaussian fits to each individual 
spectrum in the data cube, keeping only the positive- 
amplitude fits. Comparison between Gaussian fitting re¬ 
sults and conventional moment maps confirm that the 
line fluxes, line-of-sight velocities and velocity dispersions 
(cr) as resulting from Gaussian fits are all reliable and 
f ully con sistent with the maps presented in |Izumi et al.| 

y 2015aI. The difference is that our fitting scheme al- 
ws retrieving maps with a wider spatial coverage and 
recovering more extended emission, making our maps 
more suitable for the analysis of the gas kinematics. We 
further estimated the contribution of the side-lobes to 
< 15% at i? < 3", hence not large enough to distort the 


derived kinematics (see Fig 
amplitudes slightly vary for 


tne 


While the emission line 
different tracers, we note 


that their morphology remains unchanged. 

Galactocentric flux and a profiles are derived apply¬ 
ing robust statistics, ensuring reliable results with up 
to 50% tolerance level. Moreover, the above estimated 
kinematics depend on the implicit assumption that the 
data comprise a random sample from a normal distribu¬ 
tion. Noisy values or regions with stronger presence of a 
secondary component are thus accounted for by using ro¬ 
bust statistics. Accordingly, we cut each two-dimensional 
map into projected circular galactocentric rings. At each 
radius (in the plane of the galaxy) we derive the median 
value for the flux and a, followed b y the error in the de¬ 
rived m edian values as described in IHiiber fc Ronchetti I 


(20091. 

INumerous tests with varying ring widths and ring radii 
confirm that the exceptional quality and broad dynami¬ 
cal range of the ALMA data allows a stable median and 
error estimation even when median fluxes and a values 
are derived at radii a factor two smaller than the syn¬ 
thesized beam size. The smallest reliable step is thus 
0'.'25, corresponding to an increment value of 82.5 pc, 
and this will be the step size adopted throughout the 
analysis presented here. Furthermore, the median veloc¬ 
ity error across the field is 5.5 km s“^, and the median 
emission line amplitude error across the field is 0.06 mJy 
beam“^ km s“^. All these parameters remain within the 
uncertainties both for the HCN and the HCO+ lines. 

At these radial steps, we also calculate the rotation 
curve by assuming that the circular rotation is the domi¬ 
nant kinematic feature and that our measurements refer 


method describee 

in (e.g., Schoenmakers et al. 

1997 

Fathi et al. 2005 

van de Ven & h'athi |201U|) to derive the 


the line-of-sight velocity is written as the harmonic se- 
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Fig. 2.— HCN J(4 3) distribution and kinematics (top), and a representation of the noncircular velocities with a harmonic model to 

quantify their amplitudes and a final residual (data minus harmonic model) at the bottom row. The maps in the top row are results of 
single Gaussian fits to the individual spectra. The bottom left panel shows the marginal amplitudes of the non-circular motions which are 
results from the fit to the observed velocity field. The cross marks the position of the dynamical centre. The filled ellipse in the top left 
panel shows the synthesized beam size, and the dashed ellipses mark the ring-like structure at ~ 500 - 1000 pc. These maps remain 
consistent for the HCO"*" J(4 3) tracer, but slightly noisier. 

radius, fci is the best representation of the rotation curve 
(as described above) and G is the gravitational constant 
(see Fig. [^. The mass curve shown in Fi g. is in full 
agreement w ith the measur e ments listed in G^zel et al.| 


ries fco(^) + Sm=i kra{R) cos {m[il) - V'm_(7?)]), where fco 
is the systemic velocity, fci is the best fit circular velocity, 
R is the galactocentric radius, km\m>i are higher o rder 
velocity amplitudes a nd ij^rn are phase shifts (see also van 
de Ven & Fathi 2010 and Fig. Inhere). Adopting 45'^ in¬ 
clination, we hnd that the mean systemic velocity is 4860 
km s“^ and the kinematic position angle is 132°. Once 
the rotation curve is derived at each radius R, we calcu¬ 
late the angular frequency 11 = ki/R and the epieyclic 
frequency k = ^RdVfi /dR -\- curves (see Fig. 

3. RESULTS AND DISCUSSION 

Rotation dominates the kinematics of the disc and our 
data show that both the HGN and the HGO’*' lines trace 
the same dynamical entity. Henceforth, we use the HCN, 


[1 ( 1995|) and |Hicks et al. (2009), and at the ring-like 
feature, our derived mass can encompass the co mbined 
masses for the star f orming regions, as derived by |Diaz^ 


Santos et al. (2007). 


3.2. Surface Density of the Dense Gas Phase 

Considering the narrow range of the critical densities 
and the similarity of the line shapes, it is reasonable 
to conclude that the HCN and HCO^ are both emit- 
ted from the same volume (IWada et al. ||2009| |Zhang et 
S ~|[20T4t , associated with the dense cores ot star form¬ 


ing giant inolecular clouds (IKrips et al. | 2008 Greve et 


properties of the dense interstellar component, associated 

al. 120091 Andrews & Thom 

pson 120111). iVIodefing by 

with the high-mass star forming gas in this galaxy. 

iVIeijerink et af. 1 

2007) and 

Kazandjian et af. 

(2012) 


3.1. Dynamical Mass 
The low central V/a (~ 0.5) found by Hicks et al. 


[| (2009|) indicates that the random motions play a sig- 
nihcant role in the central parts of NGC 7469. Our 
Fig. [3] also confirms that the random motions have to 
be accounted for when deriving the dynamical mass. 
We assume that the HCN a rises due to symmetric 
macroscopic motions and we calculate the dynamical 
mass in the central kpc radius applying the equation 
M = R[k\ + Za^)IG (Bender et al. 1992), where R is 


lar environments and found that the densities in these 
regions vary between 10^ and 10® cm“®. In particular, 
where HCN/HCO'*' ratio is greater than unity, the most 
likely number densit y is n(H 2 ) ^ 10^ cm“®. NGC 7469 
fulfills this criterion (Izumi et al. 201^^ 


In Arp 220-like environments, Greve et af. 


found a wide range of kinetic temperatures 'l ‘k = 


(2009| 


K. Moderately star forming galaxies exhibit value s closer 
to the lower end of this temperature range (e.g., Gao fc 


Solomon [2004 Krips et al. [2008 Imanishi et al. pOlO ( 
























































































































4 


Fathi et al. 



Fig. 3. — The one-dimensional profiles for the best fit rotation 
curve ki and corresponding angular frequency Q and epicyclic fre¬ 
quency k;, and the non-circular motions up to and including the 
third order (top panel), all in the plane of the sky. Both frequency 
curves are defined in section and they are here divided by five 
for illustration purposes. The dynamical mass (thick gray curve) 
accounts for the rotation curve, the HCN a values shown in the 
bottom panel, and the assumed 45° inclination. The HCO"^ a 
profile is shown in black diamonds. The shaded area marks the 
ring-like structure. 


hence, we find it reasonable to assume = 25 K in 
NGC7469. 

Using the above density and temperature throughout 
the entire region studied here, we calculate the conversion 
factor Xhcn = 2.1 x ^ 27 Mq (K km s“^ 

pc^)“^ (Krips et al. 20081. Combined with the ratio 
HCN J(4 ^ 8) / HCN j(l ^ 0)= 0.6 ( |Viti et al. 


20141 we convert the HCN flux into molecular gas surface 
density S (in Mq pc“^). We note that the uncertainties 
in the density, temperature and HCN J(4 —?> 3) / HCN 
J(1 —> 0) are large enough to dominate the uncertainties 
in the E values illustrated in Fig. Both e S and disc 
stability parameter (as described in section [3.3[ ) depend 
highly on the assumed temperature and nurnber density. 
Decreasing the HCN J(4 —)■ 3) / HCN J(1 —>■ 0) will 
linearly increase the density values, while increasing the 
temperature will decrease the conversion factor. 

Given that the active galactic nu cleus in NGC 7469 


only influences the inner < 40 pcs (Viti et al. 


Izumi et al. ||20ihb ), our E curve is not attected by the 


2014 


corresponding radiation field. We further note that our 
measur ements are prone to be am dilution up to a factor 
5 (e.g., Sakamoto et al. TMT] ) and that this effect would 
translate into a similar level of underestimation of the E. 
Hence, we provide a rough estimate of the uncertainties 
for the E and subsequent stability parameter, using a 
factor 5 as upper value, and a factor tw o higher kinetic 
temperature (e.g., iScoville et al. ||2015|) as lower E limit 
(Fig.g). 



Fig. 4.— HCN surface density S, stability parameter Q and 
instability scale Ac (top, middle and bottom panel, respectively). 
The two dashed lines in the bottom panel mark the median value 
within each region, as indicated. The shaded area marks the ring¬ 
like structure of inte nse star formation. The uncertainties in S (as 
described in sec tion [3.2[ | are taken into account when deriving Q 
(see also section [3.3[l. 


3.3. Disc Stability Analysis 

To analyze the disc stability in the central i? ~ 1 kpc 
of NGC 7469, we assume that the velocity dispersion is 
isotropic and we derive the Toomre-like stability param¬ 
eter Q for the HCN, taking into account also the disc 


thick ness (Romeo & Wiegert 2011 


Romeo & Falstad 


20131. Accordingly, Q = ‘dKa jTWU^, where all the pa- 


rameters are described above. Romeo & Wiegert (20111 
have further shown that the dispersion relation for the 
stability parameter Q converges to Jeans dispersion rela¬ 
tion. The middle panel of Fig. [^illustrates that while the 
stability parameter Q is consistently above unity where 
no star formation is observed, this parameter decreases 
below unity at the region of the ring-like structure in 
NGC 7469. The corresponding error bars account for the 
uncertainties in the surface density E and the errors in 
the derived median tr values. 

This is a strong result, since it shows that even the 
densest gas, which is likely sunk to the disk mid-plane 
(McKee & Ostriker 20071 and least prone to pertur- 
bations, is unstable in the ring-like structure. The in¬ 
stability gives rise to the intense formation of the high- 
mass yo ung stars that have been reported in the lit¬ 
erature (jDfaz-Santos et al. 2007). Another interest¬ 
ing point here is that the HCJN y highlights the role 
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of local gravitational instabilities whi ch lead to the in- infrared emission (Gao fc Solomon~||2004 ). This correla- 


2015 

) presented a detailed discussion on the physical star-format ion 1 

Lada et al. 

2 M 2 ) 

Zhang et al. 

2014 


radii. Moreover, while gas turbulence alters the condi¬ 
tion for star-gas decoupling and increases the least stable 
wavelength, it hardly modif ies the Q parameter at scales 


larger than about 1 00 pc (Hoffmann & Romeo 2012 


Agertz et al. 2015). Given the metallicity of the cen- 
tral kpc radius of JNGC 7469, it is possible that CO-dark 
molecular gas, anisotropy of the gas velocity dispersion 
and non-axisymmetric perturbations all contribute to the 
disc stability at i? < 500 pc region, by up to 50% or 
so. However, outside this radius, the low velocity dis¬ 
persion of the HGN suggests that the contribu tion of the 
CO-dark ga s and anisotrop y may be margina l (Romeo & 
Fathi 2015). Furthermore, Sani et al. (2012) found that 
the Toomre parameter traced by HGiN is above unity at 
the centre of other Se yfert galaxies. 


The prescription of Romeo & Falstad (2013) further 
allows measuring the characteristic instability scale Ac, 
i.e. the physical scales at which a realistically thick galac¬ 
tic disc becomes locally unstable as Q drops below unity. 
It is defined as Ac = 2'ku/k and deriving the scales for 
local perturbations does not exclude the role of global 
perturbation in reshaping the disc in NGC 7469. Indeed, 
a locally stable disc can still be globally unstable to bi- 
symmetric gravitational perturbations. Moreover, the 
dynamics and evolution of disc structures depend crit- 
ically on the radial profi le of the Q stability parameter 
(|Romeo & Falstad ||2013l and references therein). 

'Lhe lower panel in Fig. shows a marked difference 
between Ac interior to and on the ring-like structure in 
NGG 7469. It is striking that in this region our measured 
^ 180 pc scale (median Ac = 0.18 kpc) is fully compara¬ 
ble to the measured sizes for the star forming complexes. 
The derived local instability scale for the densest gas is 
further supported by the observed young (5-6 Myr) and 
highly obscured {Ay ~ 10 mag) population of stars that 
coexist with an int ermediate age (15-35 M yr) less ex- 
tincted population (Dfaz-Santos et al. 2007). The young 
and massive population accounts tor ~ 1/3 of the total 
stellar mass and up to 2/3 of the infrared lumin osity in 
the ring-like structure (Dfaz-Santos et al~^|2007|. 


4. CONCLUSIONS 


Studies of star formation in discs suggest that massive 
stars are predominantly for med in the denses t cores of 
giant molecular clouds (e.g., Krumholz 2014), hence it 
is commonly assumed that HGIN and H(j(J~*~~tr ace these 
dense cores (Wu et al. 2005 Lada et al. 112012). Due to 


is likely due to the combination of a past interaction with 
the neighboring IG5283 and the presence of its bar or 
spiral arms. However, the expected non-circular motions 
are not det ected in our ALMA kinematic maps for this 
galaxy (c.f. Fathi et al. ||2005 2006 2013| ) . Our ALMA 


observations of HGIN d(4 —3) and HGU'*' J(4 —>■ 3) lines 
reveal a number of interesting clues on how local gravi¬ 
tational instabilities help forming stars in this galaxy. 

We derive the radial profiles for the disc stability pa¬ 
rameter Q and the characteristic instability scale Ac in 
NGC 7469. We find signatures of disc instability on the 
ring-like structure, inside which, the disc is stable against 
star formation. At the highly star forming ring-like struc¬ 
ture, the median local instability scale is ^ 180 pc, in 
full agreement with the sizes of individual star forming 
complexes seen in Hubble Space Telescope images. Our 
stability analysis accounts fo r a disc structure with f inite 


thickness, and as discussed in Romeo & Fathi (2015) and 


the high critical densities, a tighter correlation is found 
between the intensities of these molecules and the far 


Agertz et al. (2015), turbulence could alter the condi¬ 
tion tor star-gas decoupling and increases the least stable 
wavelength, but hardly modifies the stability parameter 
Q at scales larger than about 100 pc. 

Our robust results further show that short-wavelength 
approximation is satisfied in the region studied here and 
that the ring-like structure could in essence be only a 
tightening of the galactic-scale non-axisymmetric spiral 
arms. It will be interesting to apply this analysis to the 
warmer, less dense and more diffuse CO gas, how ever, 
as it was pointed out by |Romeo & Falstad | (|2013|) such 
a comparison must be based on identical data quality 
and resolution. At the moment we have no CO data at 
hand that match the resolution, sensitivity, and spatial 
coverage of the ALMA observations analyzed here. 
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